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ABSTRACT 

We infer the emission positions of twin kilohertz quasi-periodic oscillations (kHz 
QPOs) in neutron star low mass X-ray binaries (NS-LMXBs) based on the Alfven 
wave oscillation model (AWOM). For most sources, the emission radii of kHz QPOs 
cluster around a region of 16 — 19 km with the assumed NS radii of 15 km. Cir X-1 
has the larger emission radii of 23 ^ 38 km than those of the other sources, which 
may be ascribed to its large magnetosphere-disk radius or strong NS surface magnetic 
field. SAX J1808.4-3658 is also a particular source with the relative large emission 
radii of kHz QPOs of 20 ~ 23 km, which may be due to its large inferred NS radius 
of 18 ^ 19 km. The emission radii of kHz QPOs for all the sources are larger than the 
NS radii, and the possible explanations of which are presented. The similarity of the 
emission radii of kHz QPOs 16 — 19 km) for both the low/high luminosity Atoll/Z 
sources is found, which indicates that both sources share the similar magnetosphere- 
disk radii. 
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1 INTRODUCTION 

The launch of the Rossi X-ray Timing Explorer (RXTE) 
has led to the discovery of kilohertz quasi-periodic oscilla¬ 
tions (kHz QPOs ) in neutron star (NS) low mass X-ray bi- 
naries (LMXBs) (Ivan der Klis et ^1 19961 : IStrohmaver et al.l 
Il99d) . These high-frequency QPOs usually appear in pairs 
(upper V 2 and lower u\) at frequencies from a few hundred 
Hz to more than 1 kHz, as shown i n Atoll sources, Z sources 
(see lHasinger fc van der Klidll989l for the definition of Atoll 
and Z sources) and the accreting mil lisecond X-Ray pul¬ 
sars (AMXPs) (see Ivan der Klij l2006l for a review). It is 
suggested that kHz QPOs might be used to test the General 
Relativity in a strong gravit ational field regime l|Miller et al.l 
1 19981 : IStella fc Vietrii[l99^. and help to constrain the NS 
Mass — Radius relation ( Miller et al.lligg'sl ). 

The frequencies of the kHz QPOs are strongly cor¬ 
related with other timing and spectral features, such as 
the photon index es of the power-law component of the 
energy spectrum dKaaret et al.l Il998l ). the noise features 
(|Eord fc van der KlisI Il998d ). the positions in the X-ray 
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color-color diagra m (e.g.. IW iinands et al.l Il997bl). the X - 
ray luminosities d M4ndez et_alj _[l999b ; lFordet_^ 2000jb 
the l ow frequency QPOs ( Psaltis et al. 1 19991 : Belloni et al.l 
I 2 OO 2 I ). where the high-/low-frequency correlation is similar 
to those in black ho le candidates and white dwarf cata¬ 
clysmic variables (e.g.l^saltis et al.lll999l : iBelloni et al.l[200^ : 

I Warner fc Woudtll2002l : Im auc hel l2(in2l) . Besides, the quality 

factors and rms amplitudes of the kHz QPOs can vary as a 
function of the QP O frequencies (e.g.. iMendez et al.ll200ll : 
iBarret et al.l[^05bl ). 

There is currently no consensus as to the origin of these 
QPOs. So me interpret a tions, such as the bea t freq uency 
model by I Miller et al.l (Il998l): Lamb_&_^ill^ (I 2 OOI ) ari d 
the resonance model by iKluzniak fc Abramowicd ( 2nnlh : 


lAbramowicz et H] (I2nn3al lbh. encounter difficulties when in¬ 
terpreting t he following observationa l characteristics of the 
kHz QPOs dBelloni et al.l 200^ 2007h. W hile the relativistic 
precession model ( Stella fc Vietril 19991 : IStella et al.l Il999l ) 
fits the Au vs.1'2 relation to the kHz QPO data quite well. 
However, it requir es the NS mass of 2 solar mass (Mq) 
(|van der Klisll2n0(i ). which is relative hi gher than the actual 

measurement results dWang et aLll201^ . _ 

_ The Alfven wave oscillation model (AWOM) bv IZhan^ 

d2004h makes a good description of the kHz QPO data 
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dBelloni et al.l[200^ . The model predicts the relative emis¬ 
sion positions of the kHz QPOs {X = Rjr, the ratio between 
the NS radius and the kHz QPO emission radins), which can 
help to infer the emission radii of kHz QPOs, then further 
infer the magnetosphere-accretion disk structure and evalu¬ 
ate the possibility of the kHz QPO mechanism. 

In this paper we analyze the emission positions of kHz 
QPOs based on AWOM: In § 2, we introduce the AWOM 
and show the kHz QPO data used in the paper. In § 3 we 
infer the emission radii of kHz QPOs and analyze the results 
for the particular sources. In § 4 we make discussions and 
conclusions. 


2 THE MODEL AND DATA 


2.1 The magnetosphere-disk radius 


It is thought that kHz QPOs reflect the motion of mat¬ 
ter in orbit at some preferred radius in t he accretion 
disk around the neutron star in LMXBs (lMiller_gLalJ 


19981: Stella & Vietri 

19991: Osherovich & Titarchuk 19991: 

Lamb & Milledl200ll: 

Zhand|2004|). And in order to consider 


the interaction between the accretion flow and the magnetic 
field of NS, we introduce the magnetosphere-disk radius Vm, 
where the kinetic energy of the free-falling gas becomes com- 
parable to the magnetic en ergy of the NS magnetosphere 
JShapiro &i Teulkolskvl[l98^ 'l: 


Tui 


^rA = 


4 

M n1/7 

2GMM2 


W 


where f is a constant factor that is usually taken as 0.5 
dGhosh fc LambI Il979h . ta is the Alfven radius, is the 
magnetic moment of NS, G is the gravitational constant, 
M is the NS mass and M is the accretion rate at the inner 
disk boundary. 

We assume the NS magnetic field to be dipolar: 

B{r)^B,{-f (2) 

r 

where B^ is the NS surface dipole magnetic field strength, 
R is the NS radius and r is the radial distance refer to the 
center of the NS, and it can be seen that the magnetic field 
strength decreases with r®. Then we adopt the accreting NS 
mass M of 1.6 Mq based on its standard value of 1.4 Mq (the 
NS mass in LMXBs is averagely increased by abou t 0.2 M© 
on account of the accretion, see lZhang et al.ll2011^ and as¬ 
sume the NS radius R of 15 km, then the magnetosphere-disk 
radius can be written as: 


^7(km)(J^rA(^^: 


(3) 


2.2 AWOM 


In AWOM, the upper kHz QPO frequency 1^2 is assumed as 
the Keplerian orbital f requency i^k of the accretion flow at 
the preferred radius r (IZhanall2004l l: 


1^2 = Vk 


GM 

4j^2j,3 


(4) 


while the lower kHz QPO frequency ui is interpreted as the 
MHD Alfven wave oscillation frequency there. The model 



Figure 1. The schematic diagram of the emission radii of kHz 
QPOs for Atoll and Z sources. R is the NS radius that is assumed 
to be 15 km, and r is the radial emission radius of kHz QPOs with 
the range of 16 — 29 km and 16 — 38 km for Atoll and Z sources, 
respectively (see Tableland Table [2] for the details). 
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Figure 2. (a) The CDF curves of the position parameter X 
values. The abscissa and ordinate represent the X values and 
the CDF values respectively, and the arrows indicate the curves 
for AMXPs {X = 0.49 - 0.88), Atoll (X = 0.52 - 0.93) and 
Z (A = 0.40 — 0.92) sources, (b) Same as (a), but for the 
emission radius r of kHz QPOs (the corresponding ranges of r 
are 18.18 - 30.51 km, 16.15 - 28.27km and 16.34 - 37.58km for 
AMXPs, Atoll and Z sources respectively), and the shaded area 
shows the radius range of 16-19 km, which contains most of the 
data of Atoll sources (> 90%) and Z sources (> 70%) (see Table 
^and Table [2] for the details). 
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Table 1. The frequency ranges and the emission radii of kHz QPOs. 


Source (25) 

h'l 

(Hz) 

(Hz) 

III 

(km) 

References 

AMXP (2) 

SAX J1808.4-3658 

499 ~ 504 

685 ~ 694 

0.87 ~ 0.88 

- 

1 

XTE J1807.4-294 

106 ~ 370 

337 ~ 587 

0.49 ~ 0.83 

18.18 ~ 30.51 

2 

Atoll (15) 

4U 0614+09 

153 ~ 843 

449 ~ 1162 

0.52 ~ 0.88 

17.02 ~ 28.72 

3 

4U 1608-52 

473 ~ 867 

799 ~ 1104 

0.77 ~ 0.92 

16.32 ~ 19.49 

4 

4U 1636-53 

529 ~ 979 

823 ~ 1228 

0.81 ~ 0.93 

16.15 ~ 18.50 

5 

4U 1702-43 

722 

1055 

0.84 

17.78 

6 

4U 1705-44 

776 

1074 

0.87 

17.21 

7 

4U 1728-34 

308 ~ 894 

582 ~ 1183 

0.71 ~ 0.90 

16.69 ~ 21.02 

8 

4U 1735-44 

641 ~ 900 

982 ~ 1149 

0.82 ~ 0.91 

16.44 ~ 18.32 

9 

4U 1820-30 

764 ~ 796 

1055 ~ 1072 

0.86 ~ 0.89 

- 

10 

4U 1915-05 

224 ~ 707 

514 ~ 1055 

0.62 ~ 0.83 

18.02 ~ 24.09 

11 

Aql X-1 

795 ~ 803 

1074 ~ 1083 

0.88 

16.95 ~ 16.96 

12 

IGR J17191-2821 

681 ~ 870 

1037 ~ 1185 

0.82 ~ 0.88 

17.05 ~ 18.31 

13 

KS 1731-260 

898 

1159 ~ 1183 

0.90 ~ 0.91 

16.53 ~ 16.71 

14 

SAX J1750.8-2900 

936 

1253 

0.89 

16.88 

15 

XTE J1701-407 

745 

1150 

0.82 

18.40 

16 

XTE J2123-058 

847 ~ 871 

1102 ~ 1141 

0.89 0.90 

16.61 ~ 16.77 

17 

Z (8) 

Cir X-1 

56 ~ 226 

229 ~ 505 

0.40 0.64 

23.36 ~ 37.58 

18 

Cyg X-2 

516 

862 

0.77 

- 

19 

GX 5-1 

156 ~ 662 

478 ~ 888 

0.51 ~ 0.89 

16.90 ~ 29.69 

20 

GX 17+2 

475 ~ 830 

759 ~ 1079 

0.79 ~ 0.91 

16.47 ~ 19.05 

21 

GX 340+0 

197 ~ 565 

535 ~ 840 

0.55 ~ 0.83 

17.97 ~ 27.18 

22 

GX 349+2 

715 

985 

0.87 

17.16 

23 

Sco X-1 

532 ~ 902 

842 ~ 1143 

0.80 ~ 0.92 

16.34 ~ 18.74 

24 

XTE J1701-462 

502 ~ 651 

761 ~ 945 

0.82 ~ 0.87 

17.33 ~ 18.34 

25 


Notes: The second and third columns show the frequency ranges of the lower- and upper-kHz 
QPOs. The forth column shows the range of the inferred position parameter X{= R/r, see 
equation based on AWOM. The fifth column shows the range of the emission radii of kHz 
QPOs. 

^: The NS radii are assumed as 15km. _ 

Ivan Straaten et al.l l2005l. IWiinands et al.l 


200f 


References: 1. Ivan Straaten et al.l I20Q5I. IWiinands et al.l 12003: 2 

IZhang e t al.ll20 06bl; 3 lyai^tmaten et al.l 2000Uva^^^^^te^^tay|[200dJ^^^^^^t^| 200£ 

4. Ivan Straat^ et al.l l2003l. iBarre^^d [200^, ^ 


i et al.l 


Altamirano et al. l2008l. Wiinands et al. Il997al, iBhattacharvval I201C 

. Di Salvo et all l2003l. 

Jonker et alJl200Cl.lJonker et aLll2002al.lLin et al.ll201ll: b.lMarkwardt e1 

a al.lll999l: 7.lFord et al. 

1998a 

: 8.lDi Salvo et al. 

I 2 OOII. Ivan Straaten et al.ll2002l,|Strohmaver et al.l 19961. iMigliari et al. 

20031. 

Jonker et al.ll20nC 

. iMendez & van der Klidl 1999al: 9. IWiinands et al.lll998d. Ford et al. 

1998h 

: lO.ISmale et al.ll997l: ll.lBoirin et al.l2000l: 12.lBarret et al. 

I 2 OO 8 I: 13. Altamirano et al 

20inl; 14. IWiinands X. van der Klisll997l: 15.lKaaret et all2002l: 16. 

Strohmaver et al.l2008l: 17. 

Tomsick et ahll 19991. iHoman et ahll 19991: 18. iBoutloukos et al.ll200^ 

19. 

Wiinands et al.1 1998a: 


22. Jonker et all 20021 

IWiinands et al.l 1998a. 

Jonker et allll998l: 23. 

O’Neill et al.ll2002l: 24. 

van der Klis et al.lll997 

.Ivan der Klis et all 

199f 

.iLin et all201ll.lMendez van der Klid2000l: 

25.lHoman et al.ll2007l. 

Homan et al.ll20l0l. 

Sanna et al.ll2010l. 



predicts the following relations between the NS mass density, 
the emission position of the kHz QPOs and the twin kHz 
QPO frequencies: 

r/i = 629(Hz)A-^/V®/y 1 - yr^Jg^)2/3 (5) 


A = 



( 6 ) 


p = 4.75 X (7) 

;/2 = viX-^^'^Vl + VT^ ( 8 ) 

V = R/r (9) 

where V 2 k = 1^2/1000 Hz, A is the parameter that is related 
to the NS mass density p, M is the NS mass, Re = 77/10® cm 
is the NS radius R in unit of 10 km, X is the position pa¬ 
rameter and r is the emission radius of the kHz QPOs refer 
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Table 2. The emission radii of kHz QPOs for the sources with measured NS masses^. 


Source (3) 

Measured NS mass 

Ref 

A 

R 

r 

Risco 


(Mq) 


(oc pl/2) 

(km) 

(km) 

(km) 

Cyg X-2 

1.71 + 0.21 

1 

0.62 + 0.04 

16.4 + 1.02 

21.23 + 1.60 

15.12 + 1.86 

SAX J1808.4-3658 

< 1.4 

2 

0.43 + 0.01 

18 ~ 20ll 

20 ~ 23II 

9 ~ I 2 II 

XTE 1820-30 

-L-^y_0.07 

3 

0.65 + 0.01 

14 

16.2l°;| ~ 16.71°;® 

11 4^^'^ 


Notes: The second column shows the measured NS masses. The forth column shows the NS density parameter A(= 
see equation inferred by the twin kHz QPO frequencies in Table [J The fifth column shows the NS radii 

O rtg 

inferred by the measured NS masses and A. The sixth column shows the range of the emission radii of kHz QPOs inferred 
by the NS radii and position parameter X in Table[T] The last column shows the ISCO radii (i.e. iJisco = 3i?s = 6GM/c^). 

See Table [T] for the twin kHz QPOs and position para meter X . _ 

II: The lower limit of the NS mass i s adopted as 1 (p ee |Millej[2002|; [Zhang_et__alJ[201lf). 

Reference: l. lCasares et al.]|201Cll : 2. lElebert et al.ll20()sil : 3. IShaposhnikOT'ir'TitarchuiIl2004l . 


to the center of the NS. According to AWOM, the NS mass 
density parameter A and the position parameter X can be 
inferred by the twin kHz QPO frequencies. 


2.3 The sample of published kHz QPO 
frequencies 

We search the literature for all published instances of kHz 
QPO frequencies, and constrain the sample to the detection 
of two simultaneous kHz QPO peaks. The sample includes 
416 pairs of twin kHz QPOs from 25 NS-LMXB sources, 
where the 23 pairs are taken from the accreting millisecond 
X-Ray pulsars, and the 190/203 ones from Atoll/Z sources. 
The adopted sources with the ranges of the lower and upper 
kHz QPO frequencies and references are reported in Table[T] 
In addition, three sources in the sample have been measured 
with the NS masses, which are shown in Table [21 


3 THE EMISSION POSITIONS OF KHZ QPOS 

We infer the emission radii of kHz QPOs based on AWOM. 
For the 22 sources without the measured NS masses (see 
Table [T}, we take the following calculation steps: 

(1) For a certain pair of kHz QPO frequencies (izi and W, 
we first calculate the value of the position parameter X 
by solving equation (|S|: 

j/2 = viX-^/'^^l + VWWC. 

(2) Then we infer the corresponding emission radius r of 
the kHz QPOs with the X value by equation di): r = 
RjX, where the NS radius R is assumed to be 15 km. 

While for the three sources with the measured NS masses 
(Cyg X-2, SAX J1808.4-3658 and XTE 1820-30, see Table 
la, we first take the following steps to infer the NS radius R 
based on AWOM: 

(1) For each source, we fit the equation ([S]): 

1^1 = 629(Hz)A-2/34yi - 

to the twin kHz QPO frequencies and obtain the value 
of the NS density parameter A. 

(2) Then we infer the NS radius R of this source with the A 


value and the measured NS mass M by solving equation 


Next we take the following steps to infer the emission radii 
of kHz QPOs for the three sources: 


(1) For a certain pair of kHz QPO frequencies, we calculate 
the X value with vi and 1^2 by solving equation 

1/2 = + 

(2) Then we infer the corresponding emission radius r of 
the kHz QPOs with the X value and the inferred NS 
radius R by equation (|9l): r = RjX. 

The ranges of the inferred X and r values for the sources 
without and with the measured NS masses are shown in 
Table [T] and Table [2l respectively. There are some aspects 
should be noticed: 


(1) The emission radii of kHz QPOs of AMXPs, Atoll and 
Z sources are distributed in the range of 18 — 31 km, 
16 — 29 km and 16 — 38 km, respectively. Fig[T]shows the 
schematic diagram of the radial ranges of the emission 
positions of kHz QPOs for Atoll and Z sources. It is 
evidence that Z sources have the larger range of the 
emission radii of kHz QPOs than Atoll sources, which 
may result from source Cir X-1 with the large emission 
radii of 23.36 — 37.58 km (see also Fig|3]). 

(2) We show the cumulative distribution function (CDF) 
curves of the position parameter X and the emission 
radius r of kHz QPOs for AMXPs {X = 0.49 — 0.88, 
r = 18.18 - 30.51km), Atoll {X = 0.52 - 0.93, r = 
16.15 - 28.72 km) and Z {X = 0.40 - 0.92, r = 16.34 - 
37.58 km) sources in Fig|2] It can be seen that most of 
the emission radii of kHz QPOs of Atoll sources (> 90% 
of the data) and Z sources (> 70% of the data) cluster 
around the region of 16 — 19 km. 

(3) XTE J1807.4-294 (r = 18.18 - 30.51 km), 4U 0614+09 
(r = 17.02 - 28.72 km), GX 5-1 (r = 16.90 - 29.69 km) 
and GX 340+0 (r = 17.97 — 27.18 km) show the larger 
range of the emission radii of kHz QPOs from ~ 17 km 
to ~ 30 km (see Table [T] for the details) than the other 
sources, and their corresponding CDF curves of the 
emission radii are shown in FigO 

(4) We also analyze the innermost emission positions of 
kHz QPOs for all the sources in the sample. Fig|4]shows 
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Figure 3. The CDF curves of the emission radii of kHz QPOs for (a) XTE J1807.4-294 (r = 18.18 — 30.51km), (b) 4U 0614-1-09 
(r = 17.02 - 28.72 km), (c) GX 5-1 (r = 16.90 - 29.69km), (d) GX 340-1-0 (r = 17.97 - 27.18 km), (e) Cir X-1 (r = 23.36 - 37.58 km), 
see Table [T] for the details. 


the CDF curve of the minimal emission radii (16.15 — 
23.36 km) of all the 25 sources, from which it can be seen 
that 22 sources share the minimal emission radii around 
the range of 16 — 19 km, while the other three sources, 
i.e. Cir X-1, Cyg X-2 and SAX J1808.4-3658, show the 
larger ones of about 23, 21 and 20 km, respectively. 


4 DISCUSSIONS AND CONCLUSIONS 

We investigate the emission radii of kHz QPOs for the ac¬ 
creting millisecond X-ray pulsars. Atoll and Z sources based 
on AWOM, and find that most of the emission positions of 


kHz QPOs cluster around the region of several kilometers 
away from the NS surface. The details of the conclusions 
are discussed and summarized as follows: 

(1) The inferred emission radii of kHz QPOs of AMXPs, 
Atoll and Z sources based on AWOM are distributed in 
the range of 18 — 31km, 16 — 29 km and 16 — 38 km 
respectively (see Table [T] Table [5] and Figm for the 
details). Cir X-1 has the obviously larger emission 
radii of 23.36 — 37.58 km than the other sources. It 
is thought that kHz QPOs reflect the motion of mat¬ 
ter in orbit at some preferred radius r in the ac- 
cretion disk (iMiller et al.l 1 19981 : IStella fc Vietril 1 19991 : 
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Figure 4. The CDF curve of the minimal emission radii of kHz 
QPOs (16.15 — 23.36 km) of all the 25 sources in Tableland Table 
[2] The dashed line shows the position of the emission radius of 
19 km and the arrows indicate the innermost emission positions 
of source Cir X-1 (rmin = 23.36 km), Cyg X-2 (rmin = 21.23 km) 
and SAX J 1808.4-3658 (rmiii = 20km) (see Tableland Table [2] 
for the details). 


Osherovich &: TitarchukI 1 19991 : iLamb fc Miller! l2Q0ll : 
Zhandl2004 ). where the magnetosphere-disk radius there 


= r satisfies the equation 

We suggest that Cir X-1 may have the large NS sur¬ 
face magnetic field _Bs, which will cause the large 
magnetosphere-disk radius, i.e. the large ones of the 
emission radii of kHz QPOs. It is not clear what physical 
parameters determine the appearance of the detectable 
kHz QPOs, nor why they only occur in the certain range 
of 1 20 km away from the NS surface, and we suspect 

that it may be related to the accretion environment of 
LMXBs, which need the further study of the accretion 
state of the systems. 

(2) XTE J1807.4-294, 4U 0614-h09, GX 5-1 and GX 340-f0 
show the larger range of the emission radii of kHz QPOs 
from ~ 17km to ~ 30km (see Table [T] and Fig®!, im¬ 
plying these sources may have the more suitable phys¬ 
ical environment to produce kHz QPOs. Less bright 
source XTE J1807.4-294 (r = 18.18 ~ 30.51km) and 
4U 0614-1-09 (r = 17.02 ~ 28.72 km) have the similar 
large range of the emission radii with the more bright 
source GX 5-1 (r = 16.90 ~ 29.69 km) and GX 340-1-0 
(r = 17.97 ~ 27.18 km), which implies the emission con¬ 
dition of kHz QPOs may not sensitive to the mean ac¬ 
cretion rate. 

(3) It can be seen from see Table[T] Tableland Fig[2]that 
the emission radius r of kHz QPOs for all the sources 
are larger than the NS radius R {r > 16 km and R is 
assumed to be 15 km, i.e. r > R), which means that 
the accretion matter that is related to the kHz QPOs 
is at the position at least one kilometer away from 
the NS surface. There may be some interpretations for 
these phenomen a: One is the effect of the local strong 
magnetic spot (IZhang et al.l l2006al i which expels the 
magnetosphere-disk not to approach the NS surface, so 
the accretion matter may at the position away from the 
NS surface. Another explanation is that the accretion 


plasma drops onto the stellar surface due to the insta¬ 
bility when crossing the innermost stable circular orbit 
(ISCO) of NS, and the ISGO radius (i7isco = 6GM/c^) 
is the inner boundary of the accretion matter, so the 
emission radius of the kHz QPOs should be larger than 
7?isco. We try to test the ISCO effect scenario: Three 
sources in the sample, i.e. Gyg X-2, SAX J1808.4-3658 
and XTE 1820-30, have both the detected twin kHz 
QPO frequencies and the measured NS masses (see 
Table®, which can be used to calculate both the NS 
radii and ISGO radii with the following equations: 


1^1 = 629(Hz)A-2/34yi - 


A = 


M 1 
Mq W 


7?isco = &GMIc?. 

The results are shown in Table m i.e. for Cyg X-2, 


R ~ 16.4 km and i?isco ~ 15.12 km, for SAX J1808.4- 
3658, 77 ~ 18 - 20 km and 7?isco ~ 9 - 12 km, for XTE 
1820-30, R ~ 14.5 km and TJisco ~ 11.4 km. It can be 
seen that the ISCO radii of the three sources are all 
smaller than their inferred NS radii, so we suggest that 
the ISCO effect may not be the main reason to cause 
the emission radii of the kHz QPOs larger than the NS 
radii. 


(4) It can be seen from Table [T] Table ® and Fig® that 
most Atoll and Z sources (19/23) are inferred to share 
the similar emission radii of kHz QPOs of 16—19 km. We 
try to probe the NS magnitude field by the kHz QPOs: 
The luminosities of Atoll sources are less than those of Z 
sources in 2-3 magnitude orders in general, thus we as¬ 
sume that the accretion rates of Atoll sources are about 
2 magnitude orders lower than those of Z sources in av¬ 
erage. The magnetosphere-disk radius rm where the kHz 
QPOs occur satisfy the equation ®: 

« 7(km)( 

the similar rni(= 16 — 19km) and the different M will 
infer the NS surface magnetic field strengths of Atoll 
sources are about one magnitude order less than those 
of Z sources , which is consistent with the prediction of 
IZhang et al.l (l2007li . Source XTE J1701-462 has been ob¬ 
served both Atoll and Z behaviors llHoman et al. lioia), 
which can go fr om one state to the o ther one in timescale 
of months, and I Homan et all l|2010h suggest this source 
has a NS surface magnetic field strength should not 
change in that timescale. We argue that XTE J1701-462 
has a steady NS surface magnetic field and should not 
change in the short timescale, and ascribe the change of 
the Atoll and Z behaviors to the different magnetic field 
strengths at the different magnetosphere-disk radii: By 
assuming the NS magnetic field to be dipolar (see equa¬ 
tion ®): 

when the source has a instantaneously high accretion 
rate, the accretion disk will go closer to the NS sur¬ 
face and shows a smaller r^, where the magnetic field 
there will be higher, then the source shows the Z source 
behavior. On the contract, the system will show the 
Atoll source behavior with the lower magnetic field at 
the magnetosphere-disk radius when it has a instanta¬ 
neously low accretion rate. So we suggest that the Atoll 
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and Z states of the source may reflect the magnetic held 
information at the different magnetosphere-disk radii. 

(5) Cir X-1 (rmin = 23.36km), Cyg X-2 (rmin = 21.23 km) 
and SAX J1808.4-3658 (r^in = 20 km) show the large 
values of the minimal emission radii of kHz QPOs com¬ 
pared with the other sources in the sample (see also 
Table [U Table [5] and Fig[T]). As for Cir X-1 and Cyg 
X-2, we try to ascribe their larger ones of the minimal 
emission radii to share the stronger NS surface magnetic 
fields that make the averaged magnetosphere-disk radii 
to be larger. For SAX J1808.4-3658, its inferred NS ra¬ 
dius is slightly large (i? = 18 — 20 km with the NS mass 
of M = 1 — 1.4Mq, see Table [2|, so its emission radii 
of kHz QPOs are relative large although it has the sim¬ 
ilar emission parameter X values to the other sources 
(r = R/X). 
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